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ABSTRACT 

The origin and stability of a thin sheet of plasma in the magnetosphere of an accreting 
neutron star is investigated. First the radial extension of such a magnetospheric disc is 
explored. Then a mechanism for magnetospheric accretion is propose d, reconsidering 
the bending wave explored by Agapitou, Papaloizou & Terquem ( 1997| ), that was found 
to be stable in ideal MHD. We show that this warping becomes unstable and can reach 
high amplitudes, in a variant of P r ingle's radiation-driven model for the warping of 
AGN accretion discs ( |Pringle| ( |1996[ )). Finally we discuss how this mechanism might 
give a clue to explain the observed X-ray kHz QPO of neutron star binaries. 
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1 INTRODUCTION 

Neutron star binaries are at the centre of numerous inves- 
tigations since they are the laboratory of extreme physics, 
combining high magnetic fields and strong gravity. These 
effects dominate the interface between the star magneto- 
sphere and the accretion disc, with which the present work 
is concerned. Although many works have been dedicated 
to descriptions of this interface, a fully consistent model is 
still not available. In the pioneering work of |Ghosh &; L amb 

a strong anomalous resistivity was assumed in the disc; 
this allows magnetic field lines anchored in the neutron star 
and rotating with it to thread the disc over an extended 
transition region, while the gas remains in differential 
rotation. More recent works have rather considered a 
weaker resistivity, resulting in a sharp transition between 
the disc and the magnetosphere. 

In subsequent work |Spruit &; Ta am ( 1990j showed that an 
interchange instability, at the disc/magnetosphere interface, 
could allow the gas to penetrate from the keplerian disc 
into the magnetosphere, and they studied its subsequent 
fate. In particular they showed that, depending on the 
configuration of the magnetosphere, the gas could reach a 
radius where its vertical motion (along the field lines) could 
in turn become unstable to a warping mode, leading to 



magnetospheric accretion. However |Lepeltier &; Aly (1996) 
later found that the gas could form a stable magnetospheric 
disc, supported against gravity by the magnetic field in a 
structure very similar to the 'helmets' of stable cold gas 
suspended on magnetic field lines, commonly observed in 
the solar chromosphere. Furthermore they found that the 
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interchange mode was strongly stabilized by line tying, 
i.e. the fact that, since magnetospheric field lines are tied 
to the surface of the neutron star, they need to be bent 
somewhere between the surface and the disc to accomodate 
interchange motion in the disc; the required energy is large 
and is expected to stabilize the mode in most realistic 
magnetic configurations. 

Numerical works (Romanova et al. 2002 2003 Bessolaz 



et al. 2008| have addressed this question, but they are 
still dominated by a relatively strong dissipation. Thus, al- 
though they can address the extremely complex case where 
the rotation axis of the disc, the spin axis of the neutron 
star, and its dipolar magnetic moment are misaligned, and 
do show the possibility of magnetospheric accretion, this 
dissipation can be considered as modeling the effect of 
these instabilities, but they cannot study the instabilities 
themselves. 

In these works the emphasis was on possible instabil- 
ities of purely MHD origin, for which instability criteria 
were derived. In this case, when the modes are unstable 
their frequency is purely imaginary; they have a zero real 
frequency (Real(V OS c) = 0) in the gas frame, and would 
always be observed Doppler-shifted to rotate at the neutron 
star rotation frequency v*. 

In the present work we reconsider warp oscillations that are 
stable in the pure MHD context, but become unstable by a 
different process of radiative origin. 

We first present in section [2] the general setup and ge- 
ometry we consider, and we show how gas can enter the neu- 
tron star magnetosphere to form a dis c. In section [3] we show 
that the bending mode considered by Agapitou, Papal oTzou| 
h Terquem (1997) is stable in the configuration considered 
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of the order of 1 kHz) to the neutron star spin frequency 



Figure 1. Schematic representation of the model. The blue re- 
gion is the neutron star magnetosphere, pinched by the keplerian 
accretion disc. The magnetospheric disc inside the magnetosphere 
is supported against gravity by magnetic tension, while the kep- 
lerian disc is a usual, centrifugally supported, accretion disc. The 
arrows schematise the physical process involved in the model: as 
the gas reaches the inner edge of the keplerian disc it is picked up 
by magnetospheric field lines and slows down to the neutron star 
spin frequency v*. Supported against gravity by magnetic ten- 
sion, it is subject to the interchange instability which transports 
it inward toward the magnetic saddle point (marked by an X in 
the figure) where the concavity of the field lines changes. Radi- 
ation from the keplerian disc and the stellar surface then makes 
it unstable to a warp wave, which can trigger magnetospheric 
accretion vertically along the field lines. 



here but that radiative forcing can make it unstable, in a 



variation of the mechanism developed by |Prin gle (1996) to 
explain the warping of the external part of AGN accretion 
discs. We also discuss how this instability, since it gives the 
gas strong vertical oscillations, might provide the transition 
from radial accretion (in the keplerian disc, followed by the 
interchange mode at its inner edge) to magnetospheric ac- 
cretion along magnetic field lines, onto the neutron star. 
We note here that our approach is different from that of 
Ferreira fe Ogilvie| ([2008]) and|Ogilvie fe Dubus|p00l), who 



consider warps of the differentially rotating disc (not the 
magnetospheric disc we consider here), excited by various 
mechanisms. It is also different from the warps observed in 
numerical MHD simulations ( |Romanova et al7||2003| ) when 
the stellar magnetic field and the disc plane are misaligned. 



2 MAGNETOSPHERIC DISC 
2.1 Magnetic geometry 

Let us first describe the geometry we consider, as shown in 
figure [I] the keplerian accretion disc presses on the dipole 
magnetic field of the neutron star, so that a fraction of the 
magnetospheric field lines are pinched and have a concav- 
ity looking outward. There is necessarily a boundary layer 
Ghosh &; Lamb|1979 ) where plasma from the accretion 



&; Taam 



(typically a few times slower). As previous authors (|Spruit 
1990||Lepeltier fc Aly|1996| , and contrary to |Ghosh 



& Lamb ( 1979| , we assume that this boundary layer is very 
thin. For this and a general discussion of the physics of the 
interaction between the disc and the stellar dipole field we 



will refer to the recent works of Romanova et al. ( 2003 ) and 
Bess olaz et al.| ( |2008| , where the emphasis was on T Tauri 
stars but where the general setup is similar. 
A strong toroidal current must reside in this boundary layer, 
marking the transition from field lines anchored on the neu- 
tron star to ones anchored in the disc, and the radial Lorentz 
force from this current must, at equilibrium, balance self- 
consistently gravity and centrifugal force in the transition 
from neutron star to keplerian rotation frequency. It is im- 
portant to note here that a large fraction of the luminosity 
from the system is emitted in this boundary layer, where 
the gas is rapidly slowed down from the keplerian frequency 
at the inner edge of the disc to the neutron star rotation 
frequency, and at the surface of the neutron star. 
In the resulting geometry there must thus be within the 
magnetosphere, as shown in figure [I] a radius where the 
concavity of the field lines changes from inward to outward. 
Using cylindrical coordinates (r,<f>,z), it is easy to show 
that at this radius and in the disc midplane, the vertical 
component of the magnetic field B z starts increasing with 
radius: indeed, following a field line upward from the mid- 
plane (where B r = 0) , one sees that B r is negative on a field 
line within this radius, and positive outside this radius: thus 
dB r /dz changes sign at this radius and so does dB z /dr, 
since outside the midplane we assume to be in vacuum with 
no currents, giving VxB = 0. We call this region, marked 
by an X in figure [lj the magnetic saddle, since one easily 
shows that d 2 \B\/dz 2 is negative while d 2 \B\/dr 2 is pos- 
itive ; we call the corresponding radius r sa ddie- We do not 
expect that force- free currents, which could exist if the mag- 
netosphere is filled by a low-density plasma, would change 
this geometrical property. 



2.2 Interchange instability 



Spruit & Taam (1990); Lepeltier & Aly (1996) and Agapi- 



tou, Papaloizou &; Terquem (1997) have studied the MHD 



equilibrium and stability of the magnetospheric disc formed 
of plasma trapped in this configuration. One instability in- 
volves essentially radial motion, and the other one vertical 
motion. The first one is an interchange instability, classical 
in plasmas supported by magnetic fields against gravity: it 
interchanges flux tubes and the plasma they contain, and 
thus releases gravitational energy, if more plasma moves ra- 
dially in than out, i.e. if initially the mass per unit magnetic 
flux increases outward. ISpruit & Taam (1990) showed that 



this could allow plasma from the accretion disc to penetrate 
in the magnetosphere, and this gas forms what we call the 
magnetospheric disc, suspended against gravity mostly by 
the magnetic field; Lepelti er &; Aly| ( |1996| derived the stabil- 
ity condition of the interchange mode in the magnetospheric 
disc as (in our notations): 



disc is picked up by magnetospheric field lines, rapidly slow- 
ing down from the keplerian rotation frequency (typically 
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where E is the surface density in the disc (assumed infinitely 
thin), M is the mass of the neutron star, and A is a param- 
eter describing the geometry of the field line. 
In this equation the left-hand side is the radial derivative 
of the mass per unit of magnetic flux, as expected from the 
physics of the interchange instability. The right-hand side 
describes the stabilising effect of line tying, i.e. the fact that 
the field lines are anchored on the star: this implies that, al- 
though the instability involves motion essentially transverse 
to the field lines, some torsion must be applied to them as 
they are exchanged at the disc but not at the star surface. 
The energy this requires decreases or cancels the gravita- 
tional energy released in the interchange, stabilising it. A is 
infinite on open field lines, is large on the outermost mag- 
netospheric field lines, because they are long and can be 
twisted easily, and becomes ~1 and strongly stabilising on 
the innermost field lines. The presence of E in the left-hand 
side implies that the stabilizing effect is extremely strong if 
the surface density of the magnetospheric disc is low, i.e. the 
release of gravitational energy cannot compete with the nec- 
essary magnetic energy. 

Even without line- tying one can show that the interchange, 
near marginal stability, cannot transport the gas inward 
from the magnetic saddle, which is also the radius where 
A decreases sharply. We thus conclude that the interchange 
instability results in penetration of the gas from the kep- 
lerian disc, forming a tenuous disc where line tying is not 
yet efficient and accumulating where line tying becomes 
strong and quenches the instability. The exact radius where 
this happens cannot be farther in than the magnetic sad- 
dle, and depends very strongly on details of the keplerian 
disc/magnetosphere interaction, which could be obtained 
only from a self-consistent computation such as the ones 



of |Bessolaz et al. (2008), in 3 (rather than 2.5) dimensions 
in order to show the interchange mode. 
This (a sudden drop in density at the magnetospheric radius, 
and a local density maximum inside the magnetosphere) is 



indeed observed in the simulations of Bessolaz et al. (2008), 



and shown in their figure 4. However there the penetration 
of the gas within the magnetosphere is due to dissipation 
and can thus only roughly represent the action of the in- 
terchange mode, whose development is prevented by the as- 
sumed axisymmetry. The quenching of radial transport, re- 
sulting from line tying stabilisation of the interchange mode, 
is thus not taken into account. 



3 WARPING INSTABILITY 
3.1 Bending wave 

What happens beyond this was tentatively described by 
|Spruit &; Taam| ([i990) as resulting from the vertical MHD 
instability, which causes a warp of the magnetospheric disc 
and gives the gas some vertical motion: it can thus be ex- 
pected to trigger magnetospheric accretion, vertically along 



the field lines. However Lepeltier & Aly (1996) andfAgapi- 



tou, Papaloizou &; Terquem| ( |1997| ) refined the theory of this 



instability, taking into account rotation. We follow their sim- 
ple approach (an aligned rotator), where the disc is infinitely 
thin and the neutron star magnetic axis is aligned with its 
spin, so that the unperturbed system is axisymmetric and 



the magnetic field purely poloidal: B = (B r ,0, B z ). 
Using a WKB approximation for the radial structure of the 
bending wave, Agapitou, Papaloizou & Terquem (1997) find 
its dispersion relation: 
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where v v is the frequency of the wave, z/* is the neutron star 
spin frequency, Vk the vertical epicyclic frequency (equal to 
the keplerian rotation frequency) , B^ is the radial magnetic 
field at the disc surface, E the disc surface density and k and 
m are the radial and azimuthal wavenumbers of the pertur- 
bation ^ z ro exp[i(/cr + 7710) — 2i7vist] . 

In this equation, the right hand side terms are given by the 
vertical restoring forces applied to the disc when it is dis- 
placed from its equilibrium position: the gravitational force 
appears through the epicyclic frequency, and the Lorentz 
force through the last two terms which correspond respec- 
tively to magnetic tension and magnetic pressure. 
Furthermore, in equation [2] B^ is due to the azimuthal cur- 
rent in the disc, and this current in turn is responsible for 
the radial support of the gas, suspended in the gravitational 
field by magnetic tension: the radial equilibrium is main- 
tained by the balance of the verticaly integrated gravity, 
centrifugal and Lorentz forces: 

- 4vrVrE - 2B+B Z = (3) 

or 



where 4> is the gravitational potential of the neutron star, 
linked with the keplerian frequency, so that one finally gets 
( | Agapitou, Papaloizou &; Terquem|1997 ) 



BJ 

~eT 



47r 2 r(*4: - vl) 
B z 



(4) 



Thus equation |2] can be rewritten as 
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The surface density of the disc thus disappears from the 
second term in the right-hand side, which now depends 
only on the magnetic configuration. The third term, on 
the other hand, can in a first approach be neglected if the 
surface density of the gas is low enough. 

In this approach one gets the stability criterion, for the 
MHD warp wave, from the condition that the right-hand 
side of equation |5] be positive: 
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Since observations show that at the inner edge of the 
keplerian disc vk is typically 2 to 4 times higher than zy*, 
and is even higher in the magnetosphere, the second term 
is strongly stabilising so that the vertical mode becomes ac- 
tive only at a significant distance inward from the magnetic 
saddle (where dB z /dr = 0). This leaves a radial gap where 
the interchange can no more cause the gas to move in, and 
the vertical instability cannot yet cause it to move up along 
the field lines: the conventional MHD instabilities thus still 
leave a missing link between radial and magnetospheric ac- 
cretion. 
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Romanova et al. (2002) and Bessolaz et al. (2008) find in equation [5] 



their simulations that funnel flows can be initiated by the 
vertical pressure gradient in the gas. We note however that 
in this region of a neutron star magnetosphere (in contrast 
with the protostellar case with which these works are more 
concerned) the gas should be submitted to very rapid cool- 
ing, whereas in the absence of differential rotation, and as- 
suming that the interchange mode is stabilised by line tying, 
there is no known source of turbulence to heat it, so that the 
formation of a pressure-driven funnel flow may be difficult. 
In the present work we thus return to Spruit's original idea 
and present a model that could provide an alternative trig- 
ger for magnetospheric accretion. For this we will consider 
an additional mechanism that can make the warp unstable. 
Interestingly, this mechanism relies on the intense radiation 
field in the vicinity of an accreting neutron star. 



3.2 When the wave becomes an instability 

We now turn to the radiative mechanism that can make the 
bending wave unstable, permitting it to reach high ampli- 
tudes. For this we use a variant of the model presented by 
Pringle (1996) to explain the warp observed in the outer re- 



gion of the accretion disc of AGNs. In that case radiation 
coming from the inner region of the disc and pressing on the 
surface of the outer disc region was shown to make the warp 
wave unstable. Here on the other hand we are interested in 
an instability of the innermost (magnetospheric) region of 
the disc, and radiation can be assumed to come both from 
the surface of the neutron star and from the inner region of 
the keplerian disc, immediately outside the magnetospheric 
disc. We find that both can result in warp instability, and 
we present here a derivation only for illumination from the 
keplerian disc (i.e. from outside) since illumination from the 
stellar surface (i.e. from inside) is essentially the case stud- 
ied by Pringle - with the difference that our disc is in solid 
rather than keplerian rotation, at the stellar spin frequency 
v*. 

We assume for simplicity that the radiative flux is axisym- 
metric and comes from a transition ring at the magneto- 
spheric radius; integrating over angles we find that the lu- 
minosity dL pressing on a surface dS of the magnetospheric 
disc is, in the same polar grid as before: 
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where x is the radial distance between the surface element 
dS and the magnetospheric radius. Then the pressure force 
is 
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where c is the velocity of light and £ z the vertical displace- 
ment of the surface element. Adding this force in the equa- 
tions of motion we get the dispersion relation, modified from 
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To lowest order, the wave frequency will be close to 
v v « razy* + vk ■ 

Treating the luminosity contributions as perturbations in 
order to identify the trend of their consequences, one finds 
easily that 

• waves with negative k, i.e. leading waves are amplified 
by radiation (as in Pringle 's mechanism) 

• these waves propagate radially inward, from the keple- 
rian disc/magnetospheric disc interface toward the magnetic 
saddle 

• their radial group velocity is proportional to E, i.e. they 
propagate very slowly if the surface density of the magneto- 
spheric disc is weak. This leaves them ample time to be am- 
plified from small perturbations at the disc/magnetosphere 
interface. 



As in Pringle 's mechanism, waves with a large radial 
wavenumber have stronger linear amplification, but we ex- 
pect that, as in AGN discs, non- linear and viscous effects 



Ogilvie & Dubus (2001 ) and references therein) will se- 



lect the dominant mode. We note that in AGN the dominant 
observed mode has m = 1 and a small k (small winding an- 
gle), but there as in neutron star binaries modes of smaller 
wavelength might result in a weaker modulation and be more 
difficult to observe. 

Taking ^ l/x as the smallest possible wavelength, we note 
that the real and imaginary contributions due to the radia- 
tion field are of the same order in equation [8] This implies 
that a luminosity sufficient to induce strong amplification 
will also significantly decrease the real part of the frequency. 



3.3 The retrograde wave 

The dispersion relation, equation [8] also admits a retrograde 
wave obtained by taking the negative square root in the 
right-hand side: 
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Examination of this relation shows that this wave would 
have the same characteristics as the "direct" one we have 
discussed: leading waves with negative k are amplified and 
propagate inward with the same growth rate and group ve- 
locity. However we note that a retrograde wave would also be 
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possible in Pringle's model for AGNs, but does not seem to 
appear in observations. We believe that this might be due 
to the very low level, at these frequencies, of the thermal 
noise from which the unstable wave has to grow. We will 
thus ignore the retrograde wave in our discussion. 



4 TOWARDS A MECHANISM FOR KHZ 
QUASI-PERIODIC OSCILLATIONS ? 

If the radiation field is sufficient to make the warp wave 
strongly unstable, one should expect this to have observable 
consequences since it would result in the presence of a tilted 
surface, rotating at a frequency in the kHz range in the im- 
mediate vicinity of the most emissive regions of the star-disc 
system. Both obscuration of emissive regions and reflection 
of the luminosity should result in a strong modulation of 
the X-ray signal at that frequency, i.e. in a Quasi- Periodic 
Oscillation (QPO). In this section we explore the possibil- 
ity that the radiation-driven warp could explain one of the 
twin kilohertz QPOs observed in low-mass X-ray binaries 
that host a neutron star. We find that, although this could 
solve some questions raised by models of these QPOs, it re- 
sults in frequencies that do not correspond to the observed 
ones. This estimate is however still limited by the available 
models of the disc-magnetosphere interface. We will finally 
discuss possible properties of the disc-magnetosphere inter- 
face, not included in these models, that could produce warp 
frequencies in better agreement with the observed QPOs. 

4.1 Observations 

The frequencies v\ and V2 of the lower and upper kHz 
QPO can change by a factor of two in a given source 
but their separation Av = z/2 — z/i varies little and was 
initially observed to stay close to the frequency of the burst 
oscillations, believed to be the neutron star spin frequency 
v*. However Av is observed to decrease by a significant 
fraction when v\ and reach their highest values. Until 
recently it was generally thought that Av stayed close to 
either the burst oscillation frequency or its half, but after a 



new examination of the data Mendez & Belloni (2007) have 



shown that the link between the two quantities is not so 
direct, as we will discuss below. 



4.2 A new approach 

Numerous models have been proposed to explain the differ- 
ent types of QPO in LMXB. For the neutron star binaries 
kHz QPO, two main classes of models can be distinguished: 
those based on a beat frequency and those based on rela- 
tivistic precession motion. (iLamb et al.| (|1985|), |Miller et al 
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( 1998| ), |Stelia fe Vietrij fl998| ), jStella fe Vietrij (1999) ) 
We note here that a presupposition of both of these models, 
that z/2 corresponds to the keplerian rotation frequency at 
the inner disc edge, has remained widespread without fur- 
ther justification. 

In these interpretations of the kHz QPO, two points 
stand out. The first one is that in these models the most 'fun- 
damental' mode is the one associated with the keplerian ro- 
tation frequency at the inner disc edge, at frequency v 2 . We 



would thus expect the higher kHz QPO, to have a higher co- 
herence, and thus a higher Q factor (where Q = v/FWHM, 
and FWHM is the frequency width at half maximum) than 
the lower one at frequency v±. But the opposite is observed, 
i.e. the lower kHz QPO can have an extremely high Q, much 



higher than the other QPO (see e.g. Barret et al. (2005)). 

The second point we wish to address has been empha- 
sised by |van der Klis ( 2000 ) : although the beat frequency 



model fails to explain the variations of Az/, the proximity 
between Av and v* would indicate that somehow the gas 
disc seems to 'know' about the neutron star spin frequency. 
This would be very difficult to explain, since any direct con- 
nection (e.g. by magnetic field lines) between the disc and 
the neutron star or its magnetosphere would mechanically 
imply a constant Av. 

In the present work we consider the gas that has entered 
the magnetosphere, and thus rotates at v* (and will presum- 
ably end up accreted along magnetospheric field lines) . Mir- 
roring the previous argument, it would be difficult for this 
gas to know about the keplerian rotation frequency in the 
disc. However this gas feels the same gravity field that cor- 
responds to the keplerian rotation, although magnetic sup- 
port allows it to rotate solidly at the stellar frequency. This 
gravity field also defines the vertical epicyclic frequency of 
the magnetospheric gas, so that according to equation [9] the 
warp frequency v v should be close to vk + ^* for an m = 1 
mode. Thus assuming that the lower kHz QPO (rather than 
the higher one as in other models) appears at the rotation 
frequency at the inner disc edge, one could expect the warp 
to appear as a higher-frequency QPO with a frequency dif- 
ference close to v* . 

However additional terms appear in equation[9] and can sub- 
stantially change the warp frequency. The main unknown 
quantities that can affect the frequency difference Av are 

(i) the effect of the second term, due to the gradient of 
B z , in the right-hand side of equation [9] 

(ii) the distance between the magnetospheric radius, 
where we assume that the QPO at frequency v\ originates, 
and the radius where the frequency of the warp must be 
computed 

This should be obtained from a self-consistent computation 
of the magnetospheric disc equilibrium, in three dimensions 
and with low enough dissipation to show the interchange 
instability and its consequences on the accretion of the gas. 
This is far beyond present theoretical or numerical possi- 



bilities. Lepeltier &; Aly (1996) have studied simple exact 



models of such equilibrium configurations, in the case of 
a partially diamagnetic keplerian disc. They found a class 
of stable equilibria that we can use to estimate representa- 
tive values of the magnetic terms of the dispersion relation. 
We present in appendix detailed calculations of these terms, 
which let us calculate the resulting warp frequencies using 
equation [9] 

Fig [2^ shows the resulting evolution of Av as a function 
of v* in one of these models. The free parameters are the 
neutron star mass, spin and magnetic moment, the position 
tm and n of the inner edges of the keplerian and magneto- 
spheric disc, and finally the radius ro of the magnetospheric 
disc where the frequency of the wave is computed. We have 
chosen the units such as tm — 1, GM = 1 where G is the 
gravitational constant and M the neutron star mass. The 
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magnetospheric disc extends between n = .8 and tm — 1, 
and the frequency of the wave is computed at ro = .8, the 
location of the magnetic saddle. 

The general trend is that, as observed when comparing dif- 
ferent sources, Av decreases when v* increases, although 
the ratio between them is higher than observed, but this de- 
pends both on the radii chosen and on the choice we have 
made in the class of equilibria of |Lepeltier &; Aly| ( |1996| ). 



4.3 The disc/ magnet osphere interface 

Figure [2]d shows the same plot as Figure |2^l using a different 
model, where the saddle is closer (at ro = .88) and the disc is 
not fully diamagnetic, so that it is threaded by a fair amount 
of the neutron star magnetic flux. Av is decreased but still 
remains larger than v* . Obviously at this stage these results 
do not permit to consider the radiatively-driven warp as an 
explanation for the observed QPO. There are however two 
effects that could reduce Av to values compatible with the 
observations. 

The first effect is the role of the luminosity, in equation [9] It 
contributes both to the growth rate and to the real part of 
the frequency, in a ratio equal to kx. Let us thus consider a 
luminosity high enough to make the warp strongly unstable, 
with a growth rate of the order of the rotation frequency: its 
effect on the warp frequency, for a large-scale warp (kx ~ 1) 
will thus be to decrease Av by a similar amount, bringing 
it to a better possible agreement with the observations - in 
particular the fact that Av decreases at high luminosity. 
The second effect can be understood by considering the ra- 
dial profile of B z , as shown in figure [3] with different sets 
of parameters: B z first decreases outward (at low radius the 
dipole field of the neutron star is weakly affected by the 
disc), then has a minimum (the magnetic saddle) and in- 
creases toward the magnetospheric radius: this is due to the 
fact that in the basic model the disc is fully diamagnetic, 
rejecting the neutron star magnetic field: one thus has a 
strong current ring at the magnetospheric radius, explain- 
ing the sharp rise of B z which strongly increases the warp 
frequency, according to equation [9] When the diamagnetism 
(the parameter A of |Lepeltier fe Aly| ( [19 96 )) is reduced, a 
part of the stellar field threads the disc, reducing the ring 
current and the rise of B z - and finally Az/, as in figure [2j). 
Such a reduction of the diamagnetism is probably unreal- 
istic, since the stellar field lines rotating at v* in the disc 
would involve an extremely strong, and probably unrealis- 
tic, turbulent resistivity. However the model doesn't include 
a magnetic field originating in the disc itself, with magnetic 
field lines open to infinity rather than tied to the star. We be- 
lieve that the inclusion of this magnetic contribution would 
also reduce the ring current without requiring resistivity, 
and might result in a warp frequency in better agreement 
with the observed QPO frequency. 



5 CONCLUSION 

In this work we have explored a variant of the radiation- 



driven instability, found by Pringle (1996) to explain the 
warps observed in the outer region of AGN discs. In this 
variant radiation from the inner region of a neutron-star 
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Figure 2. The upper plot shows the dimensionless quantity 
v v (tq = O.Stm) — vk(tm)/v* = Av/v* in our model, when 
v *I v k(Tm) varies. Av is the frequency difference between the 
two kHz QPO, v* is the neutron star spin frequency and vk( t m) 
is the rotation frequency at the inner edge of the keplerian disc. 
The lower plot is the observational data as plotted by Mendez 
and Belloni (2007). The shape of the curves is similar, although 
the numerical values differ, as explained in the text. 



7 



m 4 




Figure 3. Vertical component of the magnetic field in the plane 
of the disc in function of radius. The straight, dash and dotted 
dash lines correspond to different values of the diamagnetic pa- 
rameter: respectively A = 0.3, 0.6, 1 



accretion disc destabilizes the gas that has entered the stel- 
lar magnet osphere and formed a magnet ospheric disc. The 
mechanism is similar, but can be expected to be even more 
effective as the radiating region and the destabilized one are 
in very close contact. 

The resulting warp wave can be expected to play an impor- 
tant role in the chain that permits the gas, first accreted 
in a conventional manner in the keplerian disc, to enter the 
magnetosphere by the interchange instability as found by 



Spruit &; Taam ( 1990), to form a stable magnetospheric disc 



as shown by Lepeltier &; Aly (1996), and eventually to ac- 
crete onto the neutron star by flowing along its magnetic 
field lines. More specifically, we have discussed how the warp 
could provide a crucial missing link between the interchange 
and magnetospheric accretion. 

We have finally discussed how this warp, forming a non- 
axisymmetric feature in the most strongly emitting region 
of the star-disc system, should produce a strong modulation 
of the luminosity by reflexion and obscuration of the accre- 
tion power. This could produce a Quasi-Periodic Oscillation 
in the kHz range but, at the present stage of models of the 
magnetosphere/ accretion disc interface, its frequency is not 
adequate to reproduce the observed QPO frequency. 
On the other hand the warp would be well adapted to ex- 
plain some other properties of the kHz QPOs, such as the 
variation of he frequency with the luminosity. This will be 
an incentive for further work that will aim at perfecting 
the models of the magnetosphere/accretion disc interface, 
in order to explore possibilities that could result in warp 
frequencies in better agreement with QPO observations. 
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Figure Al. The disc-magnetosphere configuration computed in 
the appendix. 



APPENDIX A: MAGNETIC CONFIGURATION 

We present here the model used for the calculation of the 
magnetic terms of the dispersion relation, equation [8] using 
the equilibria of Lepeltier & Aly (1996) which describe a 



neutron star surrounded by a keplerian disc and a magneto- 
spheric disc. They show that the magnetic field can written 

as 

B = B> + Bk + Bm 

where the right-hand side contains contributions from re- 
spectively the stellar dipole field, currents in the keplerian 
disc and in the magnetospheric disc. Introducing a param- 
eter A measuring the diamagnetism of the keplerian disc, 
which excludes totally (A = 1) or partially (A < 1) the 
magnetic field from the neutron star, they find a family of 
equilibria where these contributions can be written as 



B* = Bi((l - A)/ii,oo) 
B K - Bi(A/xi,r M ) 
B D B 3 (i,r M ) 



(Al) 
(A2) 
(A3) 



where tm is the inner radius of the keplerian disc and j 
the azimuthal current in the magnetospheric disc. Bi(/i, tm) 
and B3(j,rAf) can be derived from the azimuthal magnetic 
potentials A = fu^ with (in a-dimensional units, convenient 
since for us only the magnetic geometry is important rather 
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that the amplitude of the field) 

ef A 2a fr 2 X 2 - d 2 cos 2 6 . 2n , v \ , k A . 

fi 0, d) = — I z^z h sinOarctanX (A4) 

Tir \ d z X J 

where X = |l((l-£) 2 + ^cos^f - 1 (l - ^) } * , 
and 

poo 

fs(j,d) — rsinO I g(k)Ji(krsin6)e~ krcose dk (A5) 
Jo 

where J\ is the Bessel function of first kind and order 1 , and 

m = / d (£!!lM_ cos(fcs) )^M ds (A6 ) 

Jo s 

ow - ^IoW^ d ' (A7) 

j(n<r<r M ) J07 ^ (A8) 

and c is the speed of light. For our examples shown in figure 
[2] we have chosen \i\ — 10, jo — —10, Ti — 0.8, r# = 1., and 
respectively A = .8 and A = .265. Fig |A1| gives the shape of 
the corresponding magnetic field lines for A = .8. 



REFERENCES 

Agapitou, V., Papaloizou, J. & Terquem,C. 1997, MNRAS, 
631, 645 

Barret, D., Kluzniak, W., Olive, J. F., Paltani, S., Skinner, 

G. K., 2005, MNRAS, 1288, 1294 
Bessolaz N., Zanni C, Ferreira J., Keppens R., Bouvier J., 

2008, A&A, 478, 155 
Fe rreira, B.T. and Ogilvie, G.I., to be published in MNRAS 

( |http://arxiv.org/abs/0810.1636vT| 
Ghosh, P., & Lamb, 1979, ApJ, 234, 296 
Lamb, F. K., Shibazaki, N., Alpar, M. A., & Shaham, J. 

1985, Nature, 317, 681 
Lepeltier, T., & Aly, J.J., 1996, A&A, 645,665 
Mendez, M. and Belloni, T., 2007, MNRAS, 381, 790 
Miller, M.C, Lamb, F.K., & Psaltis, D., 1998, ApJ, 508, 

791 

Ogilvie, G.I., & Dubus, G., 2001, MNRAS, 320, 485 
Pringle, J.E., 1996, MNRAS, 281, 357 
Spruit, H.C., & Taam, R.E., 1990, A&A, 475,493 
Romanova M. M., Ustyugova G. V., Koldoba A. V., 

Lovelace R. V. E., 2002, ApJ, 578, 420 
Romanova, M.M., Ustyugova G.V., Koldoba A.V., Wick, 

J.V., Lovelace, R.V.E., 2003, ApJ, 616, L151 
Stella, L. & Vietri, M. 1998, ApJ, 492, L59 
Stella, L., Vietri, M., 1999, Phys,Rev.Let, 82,17:20 
van der Klis, M., 2000, Ann.Rev.A&A, 717, 760 
van der Klis, M., 2006, in Compact Stellar X-ray Sources, 

ed. W.H.G. Lewin and M. van der Klis (Cambridge Univ. 

Press) ( |astro-ph/0410551| ) 



This paper has been typeset from a TgX/ WT^X. file prepared 
by the author. 



